Energy of vanishing flow: mass-isospin dependence 



O 

(N 



00 
O 



Sakshi Gautam^, Aman D. Sood^, and Rajeev K. Puri^ 

^Department of Physics, Panjab University, Chandigarh - 1600 14, INDIA and 
^ SUBATECH, Laboratoire de Physique Subatomique et des Technologies Associees, 
Universite de Nantes - IN2P3/CNRS - EMN 
4 rue Alfred Kastler, F-44072 Nantes, France 



Introduction 

The investigation of the system size effects 
in various phenomena of heavy-ion colhsions 
has attracted a lot of attention. The col- 
lective transverse in-plane flow which reflects 
the competition between attractive and repul- 
sive interactions has been found to depend 
strongly on the combined mass of the system 
The energy dependence of the collective 
transverse in-plane flow has led us its disap- 
pearance at the balance energy (Ebai) A 
power law mass dependence (a A"^) of Ebai 
also has been reported Q. Earlier power law 
arameter r was supposed to be close to -1/3 
whereas recent studies showed a deviation 
from the above-mentioned power law [4| where 
T was close to -0.45. With the availability of 
high intensity radioactive beams at many fa- 
cilities, the effects of isospin degree of freedom 
in nuclear reactions can be studied in more 
details over a wide range of masses at differ- 
ent incident energies and colliding geometries. 
In the present work, we aim to study the ef- 
fect of isospin degree of freedom on the Ebai 
throughout the mass range. As reported in 
the literature, the isospin dependence of col- 
lective flow has been explained as the com- 
petition among various reaction mechanisms, 
such as nucleon-nucleon (nn) collisions, sym- 
metry energy, surface property of the colliding 
nuclei, and Coulomb force. The relative im- 
portance among these mechanisms is not yet 
clear. In the present study, we aim to shed 
light on the relative importance among the 
above-mentioned reaction mechanisms. 



The model 

In the IQMD model [5], the propagation is 
governed by the classical equations of motion: 



dH 
dp 



; Pi = 



dH 
dr., ' 



(1) 



where H stands for the Hamiltonian which is 
given by: 

i 

A 



The V^^, Vi^'\ Vf°'^, y™*, and F/^" 



are. 



respectively, the Skyrme, Yukawa, Coulomb, 
momentum dependent interactions (MDI), 
and symmetry potentials. The final form of 
the potential reads as @ 

[/™* «t4?n2[t5(pi-p2)2 + i]<5(ri-r2). (3) 

Here ^4 = 1.57 MeV and ^5 = 5 x IQ-'^MeV"'^ . 
A parameterized form of the local plus MDI 
potential is given by 

U = q(-^) + ^{-^r + 5lnMp/pof'^ + np/po- 
po po 

(4) 

The parameters a, /?, 7, 6, and e are listed in 
Ref. 161. 



Results and discussion 
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We have simulated the reactions 

24Mg-H24j^jg^ 58Cu+5«Cu, '^'^Kt+'^'^Kt , 
96Cd-f96Cd, 120Nd+120Nd, 135jJo+135Ho, 





140 




1 


1 1 


1 1 1 1 1 1 . 


leon 


120 




b/b^ 


= 0.35 - 0.45 J 


3 


100 






















> 


80 














-Coul. ftill(red.)N/Z ^ 










" • ( ♦ ) 1-0 










60 


L O ( O ) 1-4 










Coul. full (red.) 












-T|,^ =-0.S0±0.01 (-0.25±0.02) 










= -0.45±0.01 (-0.28±0.02) 










40 




1 




. 1 ... 1 














20 













:(b)( e )' Coul. full 



_ ( ^ ) Coul. red. 


1 

e 






< 






1 




" 



System Size (A) 



FIG. 1: (a) Ebai as a function of combined mass of 
system, (b) The percentage difference AEbai{%) 
as a function of combined mass of system. Solid 
(open) symbols are for N/Z = 1.0 (1.4). 



having N/Z — 1.0 and reactions ■^'^Ne+^'^Ne, 
58Cr+58Cr, 722n+722n, 962r+96Zr, 
i20Sn+i20Sn, and i^^Ba+^SBa, having 
N/Z = 1.4, respectively at semicentral impact 
parameter range 0.35 - 0.45. A soft equa- 
tion of state along with anisotropic standard 
isospin and energy dependent nucleon-nucleon 
cross section a = 0.8 aNN^'^'^'^ is being 
used. In Fig. 1(a), we display the Ef,ai as a 
function of combined mass of the system for 
the two sets of isobars. The solid and open 
circles represent the E(,az for systems with less 
and more neutron content, respectively. The 
calculated ^bai fall on the line that is a fit of 
power law nature (oc A"^), where r — -0.45 
± 0.01 and -0.50 ± 0.01 for N/Z = 1.4 and 
1.0, respectively. The different values of r for 
two curves can be attributed to the larger 
role of Coulomb force in the case of systems 
with more proton content. Our value of T1.4 
is equal/close to the value -0.45/-0.42 in Ref. 
[? ] both of which show deviation from the 
standard value ~ -1/3 where analysis was 
done for lighter mass nuclei only (< 200). 
However, for heavier systems, r increased to 
-0.45 [? ], suggesting the increasing impor- 
tance of Coulomb repulsion. This indicates 
that the difference in the Eba/ for a given pair 



of isobaric systems may be dominantly due 
to the Coulomb potential. To demonstrate 
the role of Coulomb, we have calculated the 
Efca; with Coulomb being reduced by a factor 
of 100. The results are displayed in Fig. 
1(a) with solid and open diamonds repre- 
senting systems with less and more neutron 
content, respectively. One can clearly see 
the dominance of Coulomb repulsion in both 
the mass dependence as well as in isospin 
effects. The value of ri.4 and ri.o are now, 
respectively, -0.28 ± 0.02 and -0.25 ± 0.02. 
Now with reduced Coulomb, the systems 
with more neutron content have less ^bai- 
This is because of the fact that the reduced 
Coulomb repulsion leads to higher Eha;. So 
the density achieved during the course of the 
reaction will be more due to which the impact 
of the repulsive symmetry energy will be 
more in neutron-rich systems, which in turn 
leads to less ^bai foi' neutron-rich systems 
and hence to the opposite trend for ri.4 and 
Ti.o for two different cases. In Fig. 1(b), we 
display the percentage difference AEbai{%) 
between the systems of isobaric pairs as a 
function of combined mass of system where 

AEbai{%) - ^^^TJ^ X 100. From figure, 
we see that the percentage difference between 
the two masses of a given pair is larger for 
heavier masses as compared to the lighter 
ones. However, this trend is not visible when 
we reduce the Coulomb (diamonds). 
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